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Abstract 



A chain with random hydrophobic-hydrophilic charges is studied in the 
presence of an interface separating a polar from a non polar solvent. Within 
a Gaussian variational approach in replica space, a transition is found, from a 
high temperature region, where the chain is delocalized, to a low temperature 
region, where the chain is localized at the interface. The transition tempera- 
ture diverges as the neutrality of the chain is approached. Our results are in 
agreement with an Imry-Ma type argument by Garel et al. The problem 
of replica symmetry breaking is also addressed, and the replica-symmetric 
solution is found to be unstable. 

PACS numbers: 05.40.+j, 36.20.-r, 87.15.-v 

Random heteropolymers are often studied in connection with the behavior of complex 
macromolecules such as proteins @-f|]. From a theoretical point of view, random heteropoly- 
mers are easier to study, since one can use many tools introduced in the framework of the 
statistical mechanics of disordered systems. In the case of proteins, it is believed |||| that 
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the main force driving the protein towards the folded state in physiological conditions is 
the interaction of the amino-acids with the polar solvent. In the folded state, hydrophobic 
(non-polar) residues are buried in the chain interior, whereas hydrophilic (polar or charged) 
ones are predominantly confined at the surface. 

Interesting questions arises when an interface separates a polar and a non-polar solvent. 
The latter may be thought to model the lipidic environment inside the cell membrane 0. It 
is indeed known that membrane proteins have a higher concentration of hydrophobic amino- 
acids than proteins in solution |7| . Moreover, random heteropolymers at interfaces have also 
a technological relevance, due to their effectiveness in the reinforcing of interfaces between 
two immiscible polymers @.[J. 

In this letter we study, both analytically and numerically, a very simple model of an 
ideal random hydrophilic-hydrophobic chain in presence of an interface separating a polar 
solvent (e.g. water) and a non-polar one (e.g. oil). The model has been firstly introduced 
by Garel et al [|TJ (recently, more specific model for membrane proteins has been introduced 
||10||). In their work it is shown, with Imry-Ma type arguments ||TT| and some analytical 
observations, that the quenched randomness in the sequence of hydrophobic charges causes 
a localization transition at sufficiently low temperature. The transition temperature becomes 
infinite in the case of a neutral chain, i.e. a neutral chain is always localized. The case of 
self-avoiding neutral chains at an asymmetric selective interface has been also studied with 
scaling arguments and Monte-Carlo simulations |T^JT5[| , yielding a similar scenario. 



We have employed a Gaussian variational approach for the calculation of the replicated 
partition function, in the same spirit as in references fT6Hl8f . Within this approach it 



should be possible to cope with the possible presence of many metastable states separated 



by high energy barriers [19[], via the usual Parisi ansatz for replica symmetry breaking (RSB) 
||20|| . On the other hand, the choice of a quadratic trial Hamiltonian does not provide the 
exponential decay of the monomer density which has been found both in homopolymer 



adsorption at interfaces ETJ and in Monte Carlo simulations of the model considered here 



15| . Nevertheless, we think the capability to explore correctly replica space, typical of 



the Gaussian approach, broadly makes up for this shortcoming. This is also justified 'a 
posteriori' by our results. 

Within a replica-symmetric ansatz, we obtain explicitly the temperature dependent lo- 
calization length for a neutral chain. Our result correctly reproduces the expected scaling 
behavior in the limit of high temperatures, whereas in the limit of low temperatures it shows 
a non trivial dependence on the statistical segment length of the chain, the only microscopic 
length in the problem. For a non neutral chain we obtain the expected localization tran- 
sition, which turns out to be continuous. Within a one step RSB ansatz, we show the 
replica-symmetric solution to be unstable, which is usually interpreted as a signal of a phase 
space decomposition in many energy valleys f20 |. 



The partition function of an ideal chain in d dimensions [p!q| , in the presence of an 
interface at x — —xq, is: 
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where (3 = 1/ K B T is the inverse temperature and Rj = (x^R^) (0 < % < N) are the 
positions of the monomers, ELi being the coordinates parallel to the interface. The trace 
operation is simply a multiple integral over all possible positions of the N+l monomers. The 
chain connectivity is implemented as usual with a harmonic potential, I being the typical 
length of the chain bond. The hydrophobic charge Q is assigned to the z-th monomer, 
and determines its interaction energy with the solvent. If the monomer is on the polar 
solvent side (say x > — Xq), it takes the energy —Q, whereas if it is on the non-polar 
solvent side (x < —xq), it takes the energy Q. Thus, when Q > (Q < 0), the monomer is 
hydrophilic (hydrophobic) and prefers the right side (the left side). The hydrophobic charges 
Q are identically distributed independent random Gaussian variables, with average Co an d 
variance A, P (Q) = exp — (Q — ( ) 2 / 2A 2 / V2nA^; Co = corresponds to a neutral chain. 
Any self-interaction between monomers, such as self- avoidance, is neglected in this simple 
approach. Hence, integration over Rjii can be carried out immediately and the problem 
becomes one-dimensional. 



We are interested in the quenched free energy per monomer f q = — 1/ (N/3) InZ {Q}, 
where 7 denotes the average over disorder. We introduce n replicas of the chain for a given 
disorder realization {Q} in order to employ the usual replica trick: f q = lim n _^ — l) l n - 
Thus, we deal with an effective one-dimensional homopolymer hamiltonian 

^2^2 N / n \ 2 
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The Gaussian variational approach consists in approximating TC n with a trial function 
TCt which is quadratic in the monomer coordinates Xi >a , for integer n. This provides the 
bound lnZ„ > \nZ t + (TCt — T^n)n t ^ where Z t = Tr^ Xi a y exp [— TCt] ■ For simplicity we 
consider the case of a ring polymer (xo, a = xn,o), in such a way to exploit the transla- 
tion invariance along the chain, and we introduce the corresponding Fourier coordinate k 
(x a (k) = J2f=i Xi,a exp (— iik)). In the thermodynamic limit (N ^> 1) we get: 
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where the variational parameters /i and a ab have been introduced, associated with single 
replicas and replica pairs, respectively 

Let us consider now macroscopic quantities such as the local density of a-monomers, 
Pa(x) = J2i $ ( x i,a — x), and the so called overlap of replicas a and b, q a b(x,x') = 
J2iS ( x i,a — x) 5 (x^ — x'), with a < b. Not surprisingly, the quadratic ansatz (0) yields 
a Gaussian structure for the average of both of them: 
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27T [(A + X a b) (Ao - A afe 

where Aq = JZ n dk [A -1 (k)] aa /2n and \ a b = f- n dk [A -1 (k)] ab /2ir. The 'relative' charac 



teristic length y/Xo — X a b directly measures the overlap degree of replicas a and b. Indeed, 



there is a very low probability of finding the monomers of replicas a and b at positions more 



distant than vXo — Xab. Thus, the larger the overlap length, the less the two replicas are 
interacting. In the extreme cases: X a b = 0, q a b(x,x') = ip a (x) pi, (x f ), the overlap length 
is the same as the single chain characteristic length, and replicas are not interacting at all. 
If instead X a b = Ao, q a b (x, x') = p a (x) 5 (x — x'), the overlap length vanishes, and the two 
replicas are interacting as strong as possible. 

From now on, we will consider Ao and X a b as the variational parameters^, with respect to 
which the free energy has to be stationary. The interface breaks the translation invariance, 
so that the position of the center of mass of the chain is a relevant degree of freedom. In the 
following we will consider the chain center of mass fixed in the origin, according to with 
the interface moving at the variable position — xq, which then becomes another variational 
parameter. 

Within the replica-symmetric ansatz, X a b = X V (a, b), the quenched free energy becomes: 
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where x = x {2/X ) 1/2 , v = A/A , h (A) = VA 2 + 1 - A - ln{ 1+V f+ i ), 
h! (A) = dh/dA, and e(v,x ) = J dffid3; 2 cxp ^ ( X ^ X2 ^ 2 1 sgn (xiy/i + v + a^Vl — v + xq\ 
sgn (xiyfl + v — — v + 5o) • The chain is localized if the characteristic length v'Ao 
and the chain center of mass xq remain finite in the thermodynamic limit N —>■ 00. The free 
energy has to be minimized with respect to the variational parameters Xq and Ao, which are 
single-replica quantities, and maximized with respect to v, which is a two-replicas parameter 



x The explicit dependence of (Ao,A a b) on (p,a a b) can be worked out easily in the specific cases 
we are going to consider. In the replica-symmetric case (i.e. a a b = o V (a, b)) we have A = 
- — a U l + (T + 2d / 1 ) an j \ _ \ _|_ — — — 1 ,,,,.). . Similar, but more complicated, equations hold 
for the case of one step RSB. 



Let us consider firstly the case of a neutral chain (£ = 0). The quenched free energy is 
invariant for the transformation x —>■ —x , and it can be easily seen that the minimum of 
the free energy with respect to xq is attained for xo, m in = 0. Minimization with respect to Ao, 
and maximization with respect to v, yield then Ao )m i n = I 2 / (2dy/2v max — l) , where 1/2 < 
w max < 1 is the unique solution of the equation [(1 + v) (2v — l)] 1 ^ 2 = 2f3 2 A 2 (1 — v) 3 ^ 2 jix . 
Thus, at any finite temperature the chain is localized, in agreement with heuristic arguments 

We analyze now the asymptotic behavior of our solution at high and low temperatures. 
At high temperature ((3 A <C 1) we find: 

<W ^ \ + ^ (/3A) 4 ; A 0imin ~ ^V3n (PAf 2 . (7) 

The characteristic length of the chain in the direction orthogonal to the interface is R± = 
y/Xo- We have recovered the scaling relation R± ~ (j3A)~ 2u coming from the Imry-Ma 
argument, with the exponent v = 1/2 of an ideal chain. The argument is expected to fail 
at low temperature, where it is not possible to neglect anymore the role of the microscopic 
length I. In the limit of low temperature (j3A ^> 1), we indeed find: 
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At zero temperature, the characteristic length R± — I /y 2d is of course of the same order as 
the microscopic length I. Increasing the temperature results in a power law correction with 
the non trivial exponent 4/3. 

We consider now the more general case of a non neutral chain ((q ^ 0). In this case one 
gets A 0im i n = I 2 1 (2d\/2v mafX — lj , as in the neutral case, and the following coupled equations 



for w max and x . 



[1 + v) (2v - 1) 
(1-v? 



There are two different regimes. The high temperature regime, f3A 2 / \(q\ < 1, where equa- 
tion ( |T0D does not admit a solution. The variational paramaters then read f max = 1/2, 
^o,min — °°) 5o,min = sgn (£ ) oo, and the chain is delocalized. The low temperature regime, 
13 A 2 / |Co| > 1, where a solution exists with t> max > 1/2, Ao jmm < oo, |a;o,min| < oo, and 
the chain is localized at the interface, though the center of mass is not coincident with 
the interface itself (xo,mm 7^ 0). Therefore, lowering the temperature the chain undergoes 
a localization transition at the critical temperature (5 C = \(q\ /A 2 . This is a second order 
continuous transition: if f3 \ j3 c the chain center of mass Xq moves continuously to infinity. 

The free energy is stationary at f max > 1/2, both for neutral and non neutral chains, 
and in the delocalized phase t> max = 1/2. As noticed previously, since v = A/A ^ 0, this 
means that replicas are interacting in a sensible way, even in the delocalized phase. This is 
in contrast with the Hartree approximation, which assumes instead no interaction between 
different replicas. Such approximation has been used in a new variational approach based 
on a non-stationary Green function, which results in an asymmetric ground state even in 
the case of a neutral chain IF3i . 



In order to ascertain the stability of the replica-symmetric solution, we have computed 
the quenched free energy within the so-called one step RSB ansatz. We divide the n replicas 
in n/m groups, each containing m replicas, in such a way that A a 5 = Ai if replicas a and h 
belong to the same group, and X a b = A 2 otherwise. Defining v = Ai/A and w = A 2 /A , the 
quenched free energy reads: 

if, (2nl 2 \ 1, rnl 1-m, . ., 2gT 
Pfq = -- jln [-j-y-h IB] - —h [A] + -X 0W h' [B] + 

+ 2/?Co erf ^fj + ((3 A) 2 [1 - (1 - m) e (v, x ) - me (w, x )}} , (11) 

where A = 2d\ (1 - v) /I 2 and B = 2d\ [m (l-w)+ (1 - m) (1 - v)] /I 2 . This free en- 
ergy has to be maximized with respect to v, w, m, and minimized with respect to Ao,Xo- 
Performing the change of variables (v, w) — > (s, t), with s = v — w and t = (1 — m) v + mw, 
the equation df q /ds = has always the solution s — 0, which is the replica-symmetric 



solution, whereas df q /dt | s =o = gives equation @. The stability of this solution depends 
on whether it is actually a maximum of the free energy, i.e. depends on the sign of 



d 2 f, 



ds 2 
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Notice that h" [A] < (the 'entropy' contribution) and e" (t,x ) = d e/dt > (the energy 
contribution). At high temperature entropy is dominant, d 2 f q /ds 2 \ s= o < 0, and s = is in 
fact a maximum. But at low temperature the energy contribution prevails, d 2 f q /ds 2 \ s=0 > 
and s = becomes a minimum, signalling that the replica-symmetric solution is unstable. 
Therefore, within the Gaussian variational approach, we find a glass transition to a phase 
characterized by the presence of many metastable states separated by high energy barriers, 



according to the usual interpretation of RSB ||20 | . If the transition takes place in the way 
sketched before, i.e. if it is continuous, the transition temperature is given by the equation 
d 2 f g /ds 2 \ s=0 = 0. In principle, however, it could be possible to have a discontinuous tran- 
sition if another local maximum is present which provides a better extremum for the free 
energy when the replica-symmetric solution is still stable. 

In the case of a neutral chain ((q = 0,xq = 0), we can easily compute the temperature 
(3t at which the replica-symmetric solution becomes unstable. The solution of equations 



d 2 f q /ds 2 \ s=0 = and dfjdt \ s=0 = yields v t = l/y/2, ffi = + 5/V2/A 2 , and 
A , t = (l 2 /2d) yJl + V2 for the square of the characteristic length at the transition. 

The occurrence of RSB in such a simple model, i.e. in a model without any self-interaction 
between different monomers, is higly surprising and counter-intuitive. In fact, the energy 
term —/3J2iLoCi s S n ( x i) alone does not provide any frustration. Each monomer simply 
chooses its preferred side and frustation arises from the competiton between energy and 
chain connectivity, which is an entropic effect. Thus, at sufficiently low temperature we 
would expect no frustration at all. 

In conclusion, we have studied the localization of an ideal copolymer chain at a selective 
interface. Within the Gaussian variational method and a replica-symmetric ansatz, an ex- 
plicit calculation of the relevant physical quantities has been possible. At high temperature 
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the results agree with the predictions based on Imry-Ma argument [0,0, whereas at low 



temperature a non trivial dependence on the microscopic bond length is present. Further- 
more, the replica-symmetric solution has been shown to be unstable within a one step RSB 
scheme. This would imply a breaking of the ergodicity of the system, in contrast with simple 
intuitive arguments. The Gaussian variational approach has been thought to be effective in 



describing correctly the physics of disordered frustrated systems [fl6Hl9|]. We deem there- 
fore a quite interesting issue to try to evidentiate RSB with Monte Carlo simulations, or to 
understand why the Gaussian variational method fails in predicting it, if this should be the 
case. 

We thank Flavio Seno for introducing us to this interesting subject, and Jort van Mourik, 
Jayanth Banavar and Alain Barrat for useful discussions. 
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